arXiv: 1506.08691 v2 [cs.CE] 19 Nov 2015 


Analytical reconstruction of 
isotropic turbulence spectra 
based on the Gaussian transform 

A. Wohlbrandt*\ N. Hu'’, S. Guerin% R. Ewert'’ 


° Institute of Propulsion Technology, Engine Acoustics Departmen t 
German Aerospace Center (DLR), MUller-Breslau-StrS, 10623 Berlin, Germany 
^Institute of Aerodynamics and Flow Technology, Technical Acoustics 
Getman Aerospace Center (DLR), Lilienthalplatz 7, 38108 Braunschweig, Germany 


Abstract 

The Random Particle Mesh (RPM) method used to simulate turbulence-induced broadband noise in several aeroacoustic applica¬ 
tions is extended to realise isotropic turbulence spectra. With this method turbulent fluctuations are synthesised by filtering white 
noise with a Gaussian filter kernel that in turn gives a Gaussian spectrum. The Gaussian function is smooth and its derivatives 
and integrals are again Gaussian functions. The Gaussian filter is efficient and finds wide-spread applications in stochastic signal 
processing. However in many applications Gaussian spectra do not correspond to real turbulence spectra. Thus in turbo-machines 
the VON Karman, Liepmann, and modified von Karman spectra are more realistic model spectra. In this note we analytically derive 
weighting functions to realise arbitrary isotropic solenoidal spectra using a superposition of weighted Gaussian spectra of different 
length scales. The analytic weighting functions for the von Karman , the Liepmann , and the modieied von Karman spectra are 
derived subsequently. Finally a method is proposed to discretise the problem using a limited number of Gaussian spectra. The 
effectivity of this approach is demonstrated by realising a von Karman velocity spectrum using the RPM method. 
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1. Introduction 

A stochastic noise signal of a certain spectral shape can be 
generated by convolution of a white noise signal by a filter ker¬ 
nel of an appropriate shape ID. 

One of the most common filter kernels is the Gaussian filter 
kernel that realises a Gaussian spectrum. The Gaussian filter 
is very simple and time efficient as it has benefitial characteris¬ 
tics: its derivatives and integrals are again of Gaussian shape; 
the filtering decouples in multi-dimensional space and fast filter 
methods are available, such as Purser 13 and Young & Van- 
Vliet a filters. 

But Gaussian spectra seldom represent the physics of turbu¬ 
lence. Here more elaborate spectra are needed, such as Kol¬ 
mogorov, VON Karman or Liepmann spectra. For these spectra 
the filter kernels are very complicated to use and they are fully 
coupled in space, as shown by Dieste and Gabard a. 

Siefert et al. a succeeded in realising the Kolmogorov 
spectrum using the superposition of Gaussian spectra weighted 
manually. Others have adopted this method for other kinds of 
spectra, e.g. just recently Gea-Aguilera et al. a and Kim et 
al. Q published their findings. Note that the here presented 
method has already been presented by the authors on confer¬ 
ences, but not derived in detail iiiii. 

The objective of this note is to provide a theoretical back¬ 


ground for determining the appropriate analytical weighting 
function by means of Gaussian transform M- The analytical 
weighting function is derived for the von Karman, the Liepmann 
and the modified von Karman spectra. Furthermore, an efficient 
method is proposed to discretise the weighting function with a 
limited number of Gaussian spectra. Suggestions are made to 
choose the number of filters and their length scales. As illustra¬ 
tion, the realised velocity spectrum using the Random Particle 
Mesh (RPM) method HI is compared to the analytically derived 
velocity spectrum. 


2. Method - Gaussian Transformation 

2.7. Turbulence spectra 

The most popular models for isotropic turbulence are the von 
Karman, Liepmann, and modieied von Karman models HU- 


2.1.1. The VON Karman Spectrum 

The VON Karman spectrum is commonly used to represent 
homogeneous isotropic turbulence. It satisfies the energy law 
distribution of Ic^ for the large eddies which contain most of the 
energy and reproduces the -5/3 - law in the inertial subrange. 
The energy spectrum is given by 
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with the mean turbulent velocity m, related to the turbulent in¬ 
tensity Tu and the mean flow velocity mq by - {Tu ■ uq)^, the 
integral length scale A, and the reduced wavenumber k defined 
by ^ = k* Ike, where k* - kA, and kg - 


2.1.2. The Liepmann Spectrum 

Liepmann determined turbulence longitudinal correlation co¬ 
efficients from measurements and found that they can be ap¬ 
proximated by an exponential law f(x) = exp(^). The result¬ 
ing model spectrum is given as mil El: 
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The Liepmann spectrum is comparable in shape to the von Kar- 
MAN spectrum. 


2.1.3. The MODIFIED von Karman spectrum 

According to Bechara ma the VON Karman spectrum can be 
modified to be representative over the entire wavenumber range 
including the dissipation subrange: 


Emik) = Ei;(k) exp 



(3) 
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with the Kolmogorov wavenumber kj — , where e is the 

specific dissipation rate and v is the eddy viscosity. 


2.2. Weighting function 

According to Ewert et al. HI Altering of a white noise held 
with a Gaussian Alter kernel of a specific length scale realises a 
Gaussian spectrum of the form 


Ecik) = 
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(4) 


For convenience we introduce a new spectrum eg such that 
its integral over the wavenumber range is one, i.e. 


Note that only the weighting function f(J., A) depends on the 
integral length scale A. In the following we drop A in the ex¬ 
pression of the weighting function / and write /(/, A) = /(/). 

A parameter cr is introduced to write Equation in a suit¬ 
able manner for Gaussian transform as defined by Alecu et 
al.mi- This parameter verifies the two following relationships; 
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According to Alecu et al., pik) is a zero-mean generic symmet¬ 
ric distribution, G(cr^) is the mixture function and ./L(k|cr^) is 
the zero-mean Gaussian distribution. They define the Gaussian 
transform ^ as the operator which transforms p{x) into G{cr^). 
The inverse Gaussian transform ^^'(G(cr^)) = p(x) is simply 
given by Equation Q. 

From Eq. 0 the weighting function /(/) is given as 


/ 
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2.2.7. Von Karman weighting function 

With the VON Karman spectrum given in Eq. 0 the left-hand 
side of Equation 0 becomes 
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The direct Gaussian Transform is given by Alecu et 
al. nni Eq.(4)]: 
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We are looking for a weighting function /(/, A) to realise an 
arbitrary spectrum e(k) of integral length scale A by means of a 
superposition of Gaussian spectra cq of length scales t 


e{k 
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Using Eq. 0 and 0 yields the following solution: 
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e(k. A) = I /(/, exp I - — I dZ. ( 6 ) 


where L ' is the inverse Laplace transform. Using the relation 
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where r(n) = (n - 1)! is the gamma function, we And for Eq. 0 
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and the weighting function for the von Karman spectrum is 
given by 



^ (13) 
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Figure 1: The weighting function f(l) is used to realise typi¬ 
cal turbulence spectra of length scale A using a superposition 
of Gaussian spectra of various length scales /. The analytical 
weighting functions for the von Karman (solid), the Liepmann 
(dashed), and the modified von Karman (dash-dot) spectra are 
shown. 


2.2.2. Liepmann weighting function 

With the Liepmann spectrum given in Eq. the left-hand 
side of Equation 0 becomes 


The direct Gaussian Transform is again given by Eq. ( [TOl l and 
we write 


Sip{(r)) = 


1 -1/110 
cr^ V 2cr2 ( ^2771 


■ exp 


-2s 
\ K-d j 



( 20 ) 


This can be solved analytically, with the partial fraction expan¬ 
sion M p.783]. The image function is defined by F{p) - 
H{p)IJ{p), with J{p) being a polynome about p. Eirst, we de¬ 
termine the inverse Laplace function to H(p) and l/7(p). We 
get the inverse Laplace function by applying the convolution 
theorem afterwards. 

The inverse Laplace function of 1 /J{p) is given by Eq. o 
as 


" r(17/6) 

and the inverse Laplace function of H(p) yields 
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This is of the form of the generalised Cauchy distribution shown 
in appendix of Ref ifTOl with v = 2.5 and b - We identify 
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The Gaussian transform of pf{x\v, b) is given in ifTOll as 
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So the Gaussian transform of p{x) is 
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Einally, this yields the weighting function for a Liepmann spec¬ 
trum as 


-'(//(p)) = L^'(exp(-2^))(f) 
10 else. 
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So using the convolution theorem 

L-' {H(p)IJ(p)) = L-\l/J(p)) * L-\H{p)) (25) 

^|r(T^(f-?)"^'’exp[ff(f-;l)] if f > ;! 
|0 else. 

Equation ( |20| ) simplifies to 
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This function is plotted in Eig.[T] 


(18) 


2.2.3. Modified von Karman weighting function 

With the modieied von Karman spectrum given in Eq. ([^ the 
left-hand side of Equation 0 becomes 
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and we find the weighting function for the modified von Karman 
spectrum as 
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This function is plotted in Eig. Erom the plot we see that the 
MODIFIED VON Karman and the von Karman weighting functions 
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are identical in a region of length scales above the Kolmogorov 
scale. From the modified von Karman weighting function a cut¬ 
off condition for the smalest relevant length scales can be de¬ 
rived as 


kd > 



(28) 


2.3. Weighting of velocity spectra 

Following Ref. ca, in isotropic turbulence the one¬ 
dimensional velocity spectra Enikf) are determined by the 
energy-spectrum function E{k): 


Euih) 





dk2dki. 
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As these integrals are independent of the integral length scale 
and the weighting function /(/) is independent of the wavenum¬ 
bers it is obvious that the realisation of arbitrary isotropic ve¬ 
locity spectra by superposition of Gaussian velocity spectra is 
straight forward given by combining Eq. (|^ and (|29|) to 


Eaikx > A) = J fikmuih, m- (30) 
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2.4. Two-dimensional turbulence 

The derivation has been performed in three-dimensional 
space. The application to two-dimensional turbulence is anal¬ 
ogous. As the axial velocity spectrum Eifki) is identical in 
2D and 3D space, the weighting functions f{l) apply without 
modifications to 2D turbulence, resulting in a different one¬ 
dimensional transverse velocity spectrum £ 22 (^ 1 ) and energy- 
spectrum function E{k) for 2D turbulence. 


3. Discretisation 

The main idea of the weighting function is to realise an ar¬ 
bitrary spectrum E{k) of length scale A with the superposition 
of a limited number of Gaussian spectra Ec of length scales /„ 
with 0 < m < M. For this we discretise the integral in Eq. 
to 

00 

E(k)^ J'f(l)Ec(k,l)dl (31) 

0 

M 

^J^f(UEc(k,lJAl„, (32) 

m=0 


with the spacing Al,„. 

Eor an efficient realisation we want M to be as small as possi¬ 
ble. As many orders of the wavelength k have to be covered an 
exponential distribution of the length scales seems natural. 
We propose to discretise by 

l^^hq'", ? = m = 0...M (33) 



(a) Normalised energy von Karman spectrum realised by superposition of 10 
weighted Gaussian spectra compared to the analytical solution. 



(b) Power Spectral Density (PSD) of the axial turbulence velocity fluctuations re¬ 
alised by the RPM method in comparison to the measurements provided by Cou¬ 
pland and the analytical solution for a von Karman velocity spectrum. 


Figure 2; The discretisation and application of the method to 
model turbulence spectra. 


with the minimum and maximum length scales /q and Im, re¬ 
spectively. To define the spacing A/^, the trapezoidal rule is 
applied; 


= (A/), = (AI)„='5^^. 

(34) 

Analytical parameter variations show that a realisation with 
five filters per order of A:-variation is already sufficient for to 
get a spectrum that is visually smooth. An example result is 
shown in Eig. The green curve indicates the analytical von 
Karman spectrum and the black curve is the realisation with 10 
weighted Gaussian spectra for two orders of k-variation. The 
set of blue curves show the M - \0 Gaussian spectra of length 
scales Im weighted with the analytical weighting function fKUm) 
of Eq. ( [T3] l. In comparison a realisation with only 5 weighted 
Gaussian spectra is plotted with the orange curve. The pink 
curve shows the original Gaussian spectrum of Eq. 0 of the 
integral length scale A. 

Note that the smallest used length scale f is chosen to re¬ 
solve the highest wavenumbers of interest, but the length scale 
does not have to be smaller than the Kolmogorov length scale 
Id - Eor low wavenumbers no additional Gaussian spec¬ 
tra are needed as this region is efficiently covered if the largest 
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length scale verifies Im > 4A. This is due to the fact that all 
investigated model spectra, including the Gaussian spectrum, 
follow the power law E{k) ~ in the low wavenumber region. 

As long as the ID model spectra can be assumed valid, the 
findings are problem independent. Nevertheless depending on 
the wanted accuracy another discretisation or distribution den¬ 
sity might be more appropriate and a quantification of the error 
is needed. 


4. Realisation in RPM 

Now we implement the analytical weighting functions in the 
RPM method of Ewert et al. m. The RPM method delivers un¬ 
steady fluctuating streamfunctions ij/ and the turbulent velocity 
fluctuations v by convolution of mutually uncorrelated spatio- 
temporal white-noise in a Lagrangian frame fm with the 
Gaussian filter kernel G(r) 


v = Vx(^=Vx J A G(x-x ')‘^(x, (35) 

realising a Gaussian velocity spectrum. The dimension of the 
problem is indicated by d (2 for two-dimensional turbulence 
and 3 for three-dimensional turbulence). The variance of the 
fluctuations is inferred by Ref. HI as 

_ Afxf^A^* 

R = i//i(x, t)tf/i{x, t) = -, (36) 

Po 


where po is the local mean flow density. In order to realise a 
target kinetic energy k, - with length scale A the source 
variance is defined as 


/s Z, I\ Kf 

^ ~ 3 ^ ’ 

From Eq. (|3^ the amplitude infers as 



(37) 


(38) 


5. Application 


The analytical weighting has been successfully applied by 
the authors to reproduce the measured inflow turbulence spec¬ 
trum used for the fundamental test case 1 of the AIAA bench¬ 
mark workshop [Sl with the RPM method. The energy spec¬ 
trum discretisation shown in Fig. is the one used for the 
benchmark. Turbulent fluctuations of an integral length scale of 
A = 8mm had to be resolved in the frequency range IQQHz < 
f < lOkHz- For a resolved von Karman spectrum in that range 
we needed M - 10 Gaussian spectra reaching from length 
scales lo = A/5 to Im - 4A. In Fig. 2b the resulting measured 
and synthesised turbulent spectra of the axial velocity compo¬ 
nent are compared; the results show close agreement with the 
target spectrum. 

As each filtering process, generating the Gaussian spectra of 
length scales is fully decoupled from the others, the method 
can be implemented thread-parallel. In this way the over-head 
for realising arbitrary spectra is easily handled by using addi¬ 
tional computational power resulting in negligible penalty in ef¬ 
ficiency compared to a computation realising a Gaussian spec¬ 
trum. 


6. Conlusion 

Analytical weighting functions have been derived to realise 
VON Karman, Liepmann, and modified von Karman spectra of in¬ 
tegral length scale A by integrating weighted Gaussian spectra 
of different length scales. 

A discretisation of this integral is necessary in real applica¬ 
tions. Using an exponential distribution it was shown that very 
few Gaussian realisations per order of frequency are enough for 
a smooth resulting spectrum. The discretisation is limited in the 
lower frequency range by I < 4A and the upper limit is either 
given by grid resolution or in case of the modified von Kar¬ 
man spectrum by the Kolmogorov length scale U = Inlkd- The 
method has been validated by generating a von Karman spec¬ 
trum with the RPM method by a superposition of 10 Gaussian 
realisations to resolve 2 orders of magnitude of the frequency 
range. 


Now, to realise the discretised model spectra from Equa¬ 
tions or the Gaussian energy spectra Ecik, Im) 

are weighted with f{lm)-Alm and superposed as seen in Eq. ( [32l l. 
Hence using Eq. ( [32l l, ( |38) and ( [jTl l the amplitude for each dis¬ 
crete weighted Gaussian energy spectrum Ec{k, Im) must be 
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Po 2^-dk, 
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